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Abstract 

A strict method is used to calculate the current- voltage characteristics of a double-stranded DNA. 
A more reliable model considering the electrostatic potential drop along an individual DNA molec- 
ular wire between the contacts is considered and the corresponding Green's Function is obtained 
analytically using Generating Function method, which avoids difficult numerical evaluations. The 
obtained results indicate that the electrostatic drop along the wire always increases the conductor 
beyond the threshold than without considering it, which is in agreement with recent experiments. 
The present method can also be used to calculate the current-voltage characteristics for other 
molecular wires of arbitrary length. 
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To explain experimental disputes, severa 
jlflj l to the density functional theory 11 



Recently, the progress of molecular wires Q| attracts much attention on their transport 
behaviors. A number of experimental groups have reported measurements of the current- 
voltage (I-V) characteristics of molecules, especially Deoxyribonucleic Acid (DNA) 
DNA has a special double-helix structure with complementary nucleotide base-pairs stack- 
ing onto each other, which could possibly be a candidate for one- dimensional electronic 
transport |y|. Intense experimental investigations have already been made on the transport 
properties of DNaI^. Using scanning tunneling microscope technique, Dunlap et al. jo] 
found that DNA is an insulator. Fink and schoenenberger [3] measured I-V characteristics 
of A-DNA ropes consisting of a few double-stranded DNA (dsDNA) molecules by low energy 
electron bombardment technique, and found a linear current-voltage relation. Kasumov et 
al. measured small number of DNA molecules, and observed proximity-induced super- 
conductivity. In particular, Porath et al. joj directly measured electrical transport through 
individual DNA molecules, and suggested DNA to be a semiconductor with a voltage gap. 

theoretical models ranging from the Hiikel model 
12l | have been developed. 
In the present paper, we investigate theoretically Porath's experiment and calculate 
the current between two electrodes which are connected by a DNA molecular wire. To use 
a simple model illustrating the basic physics, we follow the algorithm proposed by Mujica 
et al. Q]. The model Hamiltonian can be written as the sum of two terms, 

H = H^ + H\ (1) 

where is the non-interacting Hamiltonian for the electrodes (electron reservoirs) and DNA 
molecule, and H' represents the coupling of contacts and DNA molecule. Using extended 
Hiikel model, and H' can be expressed as Q| 

N 

H' = E^°M)(^l+E^°l")(«l+E^?l/)(/l+ (2) 

i a=l / 

N-1 

E /5 I + 1 I +h.c. 

a=l 

H' = E^ii I I +E^iv/ I N){f I +/I.C., (3) 

i f 

where the sum on | i) (| /)) runs over the state in the continuum of left (right) reservoirs. 
A DNA molecular wire consists of N sites, with one state per site, which is denoted by | a) . 



Vii ( ^Nf) is the tunneling parameter from the left (right) reservoir to the single electron 
state l(A^) in the molecule. T-matrix formalism of scattering theory gives differential 
conductance [l^ : 

9='^\G,n\' A^iEf)ABiEf), (4) 

where Gin is the (lA^) element of the matrix formula Green's Function G, it is related to the 
molecule's site 1 and A^. Ai and Ajv are semielliptical reservoir spectral densities defined 
through Newns' chemisorption theory jl^: 

, ^Jl - (^/27)2, I E/2-y \< 1, 

0, I |> 1, 

where E is measured from the center of reservoir energy band caused by chemisorption in 
the surface of the electrode. Vk {K = l,N) describe the strength of the chemisorption 
coupling between the DNA molecule and the reservoirs, and is the bandwidth of the 
reservoir. Through Lowdin's matrix partition technique [l5[, the electrode Hamiltonian can 
be replaced by a self-energy: 

N N-1 

H^jf = Si|l)(l|+E^°l«)(«l + E/5|«)(" + l| (6) 

0=1 a=l 

+Sjv I N){N I +h.c., 

where T^k {K = 1,N) are respectively the self-energy resulting from the coupling of the 
molecule to the left (right) electrode. The Green's Function is expressed as 

where z is a complex number, whose real part E is the energy of the transfer electrons. 

Mujica et al. Q| ignored the electrostatic drop along the molecule. Here we consider 
the linear voltage drop along the molecule, i.e., the electric potential between the electrodes 
varying linearly with distance. This assumption is a good approximation to the computed 
potential profile through the molecule between the two electrodes • Thus the energies of 
the sites are function of bias V: 

El = E,~qVo{a-l), a = (1, AT), (8) 
V, = V/{N-l), (9) 
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The effective Hamiltonian can be expressed in the following matrix form: 



H. 



eff 



Eb + J^i P 

(3 E,-qVo (3 ••• 

; ■•. ■•• 

• • • (3 Eb-{N- l)qVo + T.M 

where q is the average effective charge on each site. Gin can be obtained: 



(10) 



G 



IN 



(11) 



\zI-Heff\ 

where \zl — Heff\ denotes determinant. Then our task is to deduce the expression of \zl — 
Heff\. For convenience, we define: 



a 13 

(3 a-qVa (3 ■■■ 

: ■•. ■•. p 

■ ■ • (3 a-{n- 



(12) 



and from linear algebra knowledge, 
Gin 



[-1 



xN-l rN~1 



(13) 



Di^N — DiN-l^N — -D2,ArSi + D2^N-l^l^N 

where a = — Ep, i.e. the site's energy measured from the Fermi energy of the reservoirs 
without the electric field. Dpq{P = 1,2;Q = N — 1,N) is the determinant of the matrix 
obtained from A by taking rows and columns only in the range from P to Q Q| . 

We use Generating Function method to obtain Din- Setting An = Din, and An = Din, 
we find that An satisfies the following recursion relation: 



If we define 



then we obtain that 



An = [a-{n- l)qVo]An-l - P^An-2. 



F{x) = AnX" 

n=l 



(14) 



(15) 



F{x)[l -ax + 13'^x^] + {qVo)x^F'{x) = ax - /3V 



(16) 
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Since qVo is small, we can use a perturbation method, and regard qVo as a perturbation 
parameter. Up to the first order approximation, we have: 

F{x) = Fo(x) + (qVa)Fi(x), (17) 

Then we decompose into a sum of several fractions, and expand them to a power series, 
and obtain the A^. 



Ar, = Af + {qV,)A^^\ 



n+l , I — f5 . ^ri sn+l 



1 (g + y/g^ - 4/3^) " - (g - ^d' - 4/3^) 
2^ Vg^ - 4/?^ 



/02n— 1 \n _|_ \n 

= -^(n-l)(A^-A^)[l + n(l-a2 + a/.^J^)], 

o/i Ai — A2 

g , , 

a = — , (20) 



Ai 
A, 



2/5' 
a — 



/5 



/3 

In terms of our convention, A^ is the determinant -Di^at. Using similar steps, we can ob- 
tain iv-i, -D2,Ar, D2,N-i, and then obtain Giat. It is straightforward to perform similar 
calculations for high order perturbation. 

Then, we can calculate the current through the single molecular dsDNA using wide energy 
band approximation. We assume Ak{K = 1, A^) to be energy independent. For convenience, 
we consider that the two electrodes have identical Fermi energies, which are set to zero. For 
DNA molecules in equilibrium, the bases of DNA are neutral, while in transport process, 
there are charges introduced by the contacts, and charges in DNA will be redistributed. 
Since the average charge q on each site is small, the derivation of the actual potential drop 
along DNA molecule from the assumed linear drop can be ignored. 

In previous efforts, many papers ignored the situation that the sites' energies will vary 
along the molecular wire because of the voltage drop in high-intensity electric field. Some 
papers did consider this effect (e.g., Ref. but they did not obtain an analytic expres- 

sion of the differential conductance. The present paper achieves this task. In general, the 



modified Hamiltonian is more complex. Since the Green's Function must be obtained by 
calculating the inverse matrix, it is not a trivial task, especially when the matrix is large. 

Using the linear voltage drop approximation, we first obtain the explicit expression of Giat 
up to the first order perturbation, and get a more reasonable current-voltage relation. We 
show the I-V curve in Fig.l for a 4-nucleotide DNA, and Fig. 2 for a 30-nucleotide DNA. We 
find when considering the voltage drop in the molecule, the current-voltage curves always 
become much steeper near the threshold. This is in accordance with the experiment of 

n 

Porath et al. [9]: some I-V curves in their paper are much steeper near the threshold. Our 
result is intuitional, because voltage drop effect makes the electron transport more easily 
in strong electric field. Moreover, Fig. 3a of Porath's paper contains two different I-V 
characteristics of the same 30-base pair DNA sample. The difference between these two 
curves may be the result of a sudden change in the DNA wire, possibly a conformational 
change. In transport process, the charge distribution along DNA may depend on DNA's 
chemical nature, for example, the molecule's conformation. Although there was a similar 
voltage drop in the sample of both cases, the observed electrostatic drop effect are different. 
Finally, our work can be used to calculate other molecular wires. Since the computational 
complexity of our analytical deduction is not sensitive to the number of sites, so we can 
calculate long chains. 

An issue of interest is to investigate theoretically the transport property of DNA molecules 
made of heterogeneous base sequences. We are now working on this possibility. 

In summary, we consider the electrostatic potential across the DNA, which may change 
the site energy along the DNA with the applied bias. This made the computation more 
difficult. We have proposed an analytical computational method to deal with this difficulty. 
It can be used to calculate the I-V characteristic of homogeneous DNA wires of arbitrary 
length. We can use the potential drop effect to explain recent single molecular DNA exper- 
iments. 

The authors are grateful to Guangping Gao for helpful discussions and to Haijun Zhou 
for reading earlier versions of the manuscript. We finally think Prof. J. Kong for modifying 
the manuscript. 
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FIG. 1: Theoretically predicted current- voltage characteristics of single molecular double stranded 
DNA. Parameters: /3 = -2AeV, Vi = Vn = 0.5eV, 7 = lO.OeF, = 4, and Eb is l.OeV below 
Ep. The full line is calculated with q = 0, i.e., no site's energy shift effect. The dashed curve is 
calculated with q = 0.05e. The dotted curve is calculated with q = O.le. 
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FIG. 2: Theoretically predicted current- voltage characteristic of single molecular double stranded 
DNA. Parameters: /? = -O.leV, Vi = Vn = S.OeV, 7 = lO.OeV, N = 30, and Eb is l.OeV below 
Ep- The full line is calculated with q = 0, ie., no site's energy shift effect. The dashed curve is 
calculated with q = O.Ole. The dotted curve is calculated with q = 0.05e . 
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